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ABSTRACT 



We estimate the density laws of the Galactic stellar populations as a function 
of absolute magnitude in a near-polar Galactic field. The density laws are deter- 
mined by the direct fit to photometric parallaxes from Vega photometry in the ELAIS 
(V-i ■ (a = 16''10™00",(5 = -|-54°30'00"; I = 84°. 27, b = +44°. 90; 6.571 deg^; epoch 2000) 

\ field both independently for each population and simultaneously for all stellar pop- 

\^ I ulations. Stars have been separated into different populations based on their spatial 

location. The thick disc and halo best fit by an exponential. However, the thin disc 
best fits by using a sech^ law for stars at faint absolute magnitudes, 10 < M{g ) < 11, 
I 11 < M{g ) < 12 and 12 < M{g ) < 13, whereas an exponential law for stars at 

relatively bright absolute magnitudes, 5 < M{g ) < 6, 6 < M{g ) < 7, 7 < M{g ) < 8, 

8 < M{g') < 9 and 9 < M{g') < 10. The scaleheights for the sech^ density laws 
are the equivalent exponential scaleheights. Galactic model parameters are absolute 

O I magnitude dependent: The scaleheight for thin disc decreases monotonically from 

stars at bright absolute magnitudes [M{g ) = 5] to stars at faint absolute mag- 
I nitudes [M{g ) = 13] in the range 363-163 pc, except the minimum H=211 pc at 

9 < M{g ) < 10 where sech density law fits better. Its local density is flat at bright 
absolute magnitudes but it increases at faint absolute magnitudes. For thick disc, the 

. scaleheight is flat within the uncertainties. The local space density of thick disc relative 

' to the local space density for the thin disc is almost flat at absolute magnitude inter- 

vals 5 < M{g ) < 6 and 6 < M{g ) < 7, 7.59 and 7.41 per cent respectively, whereas 
it decreases down to 3.31 per cent at absolute magnitude interval 7 < M{g ) < 8. The 
axial ratio for the halo is k =0.60, 0.73 and 0.78 for the absolute magnitude intervals 
4 < M{g ) < 5, 5 < M{g ) < 6 and 6 < M{g ) < 7 respectively, and its local space 
density relative to the local space density for the thin disc is 0.06 and 0.04 per cent 
for the intervals 5 < M{g ) < 6, and 6 < M{g ) < 7 respectively (the local space den- 
sity relative to the thin disc could not be derived for the absolute magnitude interval 
4 < M{g ) < 5 due to the lack of the local space density for thin disc for this interval). 
The simultaneous fit of all three stellar populations agrees within uncertainties with 
the most recent values in the literature. Also, each parameter is close to one of the 
corresponding parameters estimated for different absolute magnitude intervals in this 
work with one exception however; i.e. the scaleheight for thick disc is relatively small 
and its error is rather large {H = IGOtH pc). 

Key words: Galaxy: stellar content - technique: photometric - survey - methods: 
data analysis 
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Table 1. Previous Galactic models. Symbols: TN denotes the thin disc, TK denotes the thick disc, S denotes the spheroid (halo), iJe is 
the effective radius and k is the axes ratio. The figures in the parentheses for Siegel et al. (2002) are the corrected values for binajrism. 
The asterisk denotes the power-law index replacing R^. 



II (TN) 


h (TN) 


n (TK) 


H (TK) 


h (TK) 


n (S) 


Re (S) 




Reference 




(kpc) 




(kpc) 


(kpc) 




(kpc) 






310-326 
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0.76 


2.8 
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1 INTRODUCTION 

Galactic models have a long history. Bahcall & Soneira 
(1980) fitted their observations with a double component 
Galactic model, namely disc and halo, whereas Gilmore & 
Reid (1983) could succeed to fit their observations with a 
Galactic model only by introducing a third component, i.e. 
thick disc. It should be noted that the third component 
was a rediscovery of the "Intermediate Population 11" first 
described in the Vatican Proceedings review of O'Connel 
(1958). The new model is discussed by Gilmore & Wyse 
(1985) and Wyse & Gilmore (1986). Galactic models have 
been an attractive topic for many research centers, due to 
their importance: Galactic models can be used as a tool to 
reveal the formation and evolution of the Galaxy. For some 
years there has been a conflict among the researchers about 
the history of our Galaxy. The pioneering work was the one 
of Eggen, Lynden-Bell & Sandage (1962) who argued that 
the Galaxy collapsed in a free-fall time (~ 2 x 10^ yr) . Now, 
we know that the Galaxy collapsed over many Gyr (e.g. 
Yoshu & Saio 1979; Norris, Bessel & Pickles 1985; Norris 
1986; Sandage & Fouts 1987; Carney, Latham & Laird 1990; 
Norris & Ryan 1991; Beers & Sommer-Larsen 1995) and at 
least some of its components are formed from the merger or 
accretion of numerous fragments, such as dwarf-type galax- 
ies (cf. Searle & Zinn 1978, Freeman & Bland-Hawthorn 
2002, and references therein). 

The researchers use different methods to determine the 
parameters. Table 1 summaries the results of these works. 
One can see that there is an evolution for the numerical val- 
ues of model parameters. The local space density and the 
scaleheight of the thick disc can be given as an example. 
The evaluations of the thick disc have steadily moved to- 
ward shorter scaleheights (from 1.45 to 0.65 kpc, Gilmore & 
Reid 1983; Chen et al. 2001) and higher local densities (2-10 
per cent). In many studies the range of values for the param- 
eters is large. For example, Chen et al. (2001) and Siegel et 
al. (2002) give 6.5-13 and 6-10 per cent, respectively, for the 



local space density for the thick disc. However, one expected 
the most evolved numerical values for these recent works. 
That is, either the range for this parameter should be small 
or a single value with a small error should be given for it. It 
seems that they could not choose the most appropriate pro- 
cedures in this topic. In fact, we cited in our previous paper 
(Karaali, Bilir & Hamzaoglu 2004, hereafter KBH) that the 
Galactic model parameters are mass dependent. Absolute 
magnitude is reasonable proxy for mass, therefore they vary 
at different absolute magnitude intervals. Hence, the param- 
eters cited by the researchers up to recent years which axe 
based on star counts cover the range of a series of parame- 
ters corresponding to different absolute magnitude intervals, 
therefore either their range or their errors axe large. Addi- 
tionally, as it was cited in our previous paper (KBH), sech^ 
density law fits better to the observed density functions for 
stars with absolutely faint magnitudes, 10 < M{g ) < 13, 
for the thin disc. We aim to use these experiences in the in- 
vestigation of this field and compare the results with those 
obtained in the field SA 114, almost symmetric relative to 
the Galactic plane. It should be noted that evaluations of 
photometric parallax (Gilmore & Reid 1983, Reid & Ma- 
jewski 1993, Siegel et al. 2002) have usually broken the fits 
down by absolute magnitude ranges. More importantly, the 
Besangon group (e.g. Robin et al. 1996) uses very sophis- 
ticated models that create multiple thin disc populations 
through population synthesis. This is a much more elegant 
and nuanced way of fitting star count parameters. However, 
we should mention that the method of photometric parallax 
is, by necessity, a simplified way of evaluating star counts. 



2 THE DENSITY LAW FORMS 

Disc structures are usually parameterized in cylindrical co- 
ordinates by radial and vertical exponentials. 
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Di{x,z) =niexp{—z/Hi)exp{-{x — Ro)/hi) (1) 

where z is the distance from Galactic plane, x is the planar 
distance from the Galactic center, Ro is the solar distance 
to the Galactic center (8.6 kpc, Buscr ct al. 1998), Hi and 
hi are the scaleheight and scalelength respectively, and rn is 
the normalized local density. The suffix i takes the values 1 
and 2, as long as the thin and thick discs are considered. It 
should be noted that the sophisticated models of Besangon 
and others use multiple thin discs to account for the range 
of populations (e.g. Robin et al. 1996). A similar form uses 
the sech^ (or sech) function to parameterize the vertical dis- 
tribution for the thin disc, 

Di(x,z) = nisech^{—z/Hi)exp{—{x — Ro)/hi). (2) 

As the secans hyperbolicus is the sum of two exponentials. 
Hi is not really a scaleheight but has to be compared to Hi 
by dividing it with 2: Hi — Hi/2. We would like to men- 
tion that the reason of using a scch'^ law is due to theoretic 
analysis which indicate that the density laws should follow 
a sech^ law for an isothermal sheet. 

The density law for the spheroid component is parame- 
terized in different forms. The most common is the dc Vau- 
couleurs (1948) spheroid used to describe the surface bright- 
ness profile of elliptical galaxies. This law has been depro- 
jected into three dimensions by Young (1976) as 



3 THE PROCEDURE USED IN THIS WORK 

In this work, we used the same procedure cited in our previ- 
ous paper (KBH), i.e. we compared the derived and theoret- 
ical space densities per absolute magnitude interval, in the 
vertical direction of the Galactic plane for a large absolute 
magnitude interval 4 < M{g ) < 13, down to the hmiting 
magnitude g^■, — 20.5: (i) we separated the stars into different 
populations by their spatial position, as a function of both 
absolute and apparent magnitude; (ii) we tried the exponen- 
tial and scch^ laws for comparison of the derived and the- 
oretical space densities for the thin disc and we found that 
a sech^ law worked better at magnitudes 10 < M{g ) < 13 
whereas an exponential density law favors at magnitudes 

4 < M{g ) < 10. This was also the case in our previous pa- 
per (KBH); (iii) we derived model parameters for each popu- 
lation individually and for each absolute magnitude interval 
we observed their differences; and (iv) the model parameters 
were estimated by comparison of the derived vertical space 
densities with the combined density laws (eqs. 7 and 9) for 
stars of all populations. In the last process, we obtained two 
sets of parameters: one for the absolute magnitude interval 
4 < M{g') < 10 and the other 4 < M{g ) < 13. As we 
argued in our previous paper, the different behavior of the 
faint stars may produce different values and large ranges for 
parameters derived in starcount studies. 



DsiR) = ris exp[-7.669{R/Rey^^]/(R/Rey^\ (3) 

where R is the (uncorrected) Galactocentric distance in 
spherical coordinates. Re is the effective radius and Us is 
the normalized local density. R has to be corrected for the 
axial ratio k = c/a, 

R=[x' + {z/Kf]'/\ (4) 
where, 

z = r sin b, (5) 

X = [Ra + cos^ b — 2Ror cosbcosl]^^^ , (6) 

r being the distance along the line of sight and, 6 and / the 
Galactic latitude and longitude respectively, for the field un- 
der investigation. The form used by the Basle group is inde- 
pendent of effective radius but is dependent on the distance 
from the Sun to the Galactic centre: 

DsiR) = ris ea;p[10.093(l - R/ Rof'^]/{R/ Ro)"^; (7) 
and alternative formulation is the power law, 

D,{R)^nJ{aZ + R'^) (8) 

where Oo is the core radius. 

Equations (1) and (2) can be replaced by eqs (9) and 
(10) respectively, as long as the vertical direction is consid- 
ered, where 

Di{z) = niexp{—z/Hi), (9) 
Di{z) = nisech^{-z/H'i). (10) 



4 THE DATA AND REDUCTIONS 
4.1 Observations 

The ELAIS field {a = 16''10'"00", <5 = +54°30'00"; I = 
84°. 27, b = -|-44°.90; 6.571 deg^; epoch 2000) was measured 
by the Isaac Newton Telescope (INT) Wide Field Camera 
(WFC) mounted at the prime focus (//3) of the 2.5-m INT 
on La Palma, Canary Islands, during seven observing runs, 
namely 1999 April 17; 1999 June 7, 9; 1999 July 16-22; 1999 
August 1-3, 7, 10, 13, 16-17; 1999 September 7-9; 1999 Octo- 
ber 3-5, and 2000 June 24-25. The WFC consists of 4 EVV42 
CCDs, each containing 2k x 4k pixels. They are fitted in a L- 
shaped pattern which makes the camera have 6k x 6k pixels, 
minus a corner of 2k x 2k pixels. The WFC has 13.5/i pixels 
corresponding to 0.33 arcsec pixel"^ at the INT prime focus, 
and each covers an axea of 22.8 x 11.4 arcmin^ on the sky. 
This field contains 54 sub-fields and each sub-field covers 4 
CCDs with a total area of 0.29 deg^. Therefore, the total 
area of each telescope pointing is 54 x 0.29 deg^ minus the 
overlapping area. In our work, the data of only 33 sub-fields 
could be used. Hence, the area of the field investigated is 33 
X 0.29 deg minus the overlapping area = 6.571 deg^ With 
a typical seeing of 1.0-1.3 arcsec on the INT, point objects 
are well sampled, which allows accurate photometry. 

Observations were taken in five bands (m/jgo; 5 > ^ > * i 
ZfiQo, where RGO denotes the Royal Greenwich Observa- 
tory) with a single exposure of 600 s to nominal 5a limiting 
magnitudes of 23, 25, 24, 23, and 22 respectively (McMahon 
et al. 2001). However, the limiting magnitudes are brighter 
when stars are considered only. In our work, we determined 
the Qq limiting magnitude for stars by estimating from the 
star count roll over in Fig. 2 as go = 20.5. Magnitudes are 
put on a standard scale using observations of ELAIS stan- 
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dard system corresponding to Landolt system ^, taken on 
the same night. The accuracy of the prehminary photomet- 
ric cahbrations is ±0.1 mag. The CCD observations are re- 
duced to the magnitudes by the INT WAS group. Three sets 
of two-colour diagrams, i.e. (m —g,g —r),{g —r,r —i) 
and (r — i ,i — 2 ), for 21 sub- fields show considerable de- 
viations due to bad reduction hence, we left them out of the 
program. The following processes have boon applied to tho 
data for the remaining 33 sub-fields to obtain a sample of 
stars with new data available for a model parameterization. 

4.2 The overlapping sources, de-reddening of the 
magnitudes, bright stars, and extra-galactic 
objects 

The data of ELAIS field are provided from tho Cambridge 
Astronomical Survey Unit (CASU) ^. In total, there are 
17041 sources in 33 sub-fields in the ELAIS field. It turned 
out that 3027 of these sources aro overlapped, i.e. their an- 
gular distances are less than 1 arcsoc to any other source. 
We omitted them, and so the sample reduced to 14014. The 
E{B — V) colour excess for the sample sources are evaluated 
by the procedure of Schlegel, Finkbeiner & Davis (1998). 

The E{B — V) colour-excess contours for the field are 
given in Fig. 1 as a function of Galactic latitude and longi- 
tude. Then, the total absorption Ay is evaluated by means 
of the well known equation 



For Vega bands we used the R\/Rv data of Cox (2000) 
for the interpolation, where A = 3581, 4846, 6240, 7743, 
and 8763A, and derived R\ from their combination of this 
with Ay (see Table 2, KBH). Finally, the dereddened Mq, go, 
To, «oi and Zq magnitudes were obtained from the original 
magnitudes and the corresponding R\. 

The histogram for the dereddened apparent magnitude 
5o (Fig. 2) and the colour-apparent magnitude diagram (Fig. 

^ htt p : // www . ast . cam. ac. uk / ~ wf csur/ technical / photom / 
^ http://www.ast.cam.ac.Uk/~ wfcsur/release/elaiswfs/ 




-0.S 0.5 1 1.5 2 

Figure 3. Colour-appajrent magnitude diagram for the original 
sample. The shaded areas for the regions that correspond to each 
stellar population is indicated. 

3) shows that there is a large number of saturated sources 
in our sample. Hence, we excluded sources brighter than 
go = 17 (this bright limit of apparent magnitude matches 
with tho one claimed in our previous paper, KBH). How- 
ever, the two-colour diagrams (w — g )o — {g — t )o and 
(5 )o — (f — * )o in Fig. 4 indicate that there are also some 
extragalactic objects, where most of them lie towards the 
blue as claimed by Chen et al. (2001) and Siegel et al. (2002). 
As claimed in our paper cited above, the star/extragalactic 
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object separation based on the "stellarity parameter" as re- 
turned from the SExtractor routmcs (Bertin & Arnouts 
1996) could not be sufScient. We adopted the simulation of 
Fan (1999) in addition to the work cited above, to remove 
the extragalactic objects in our field. Thus we rejected the 
sources with {u — g )q < —0.10 and those which lie out- 
side of the band concentrated by most of the sources. Af- 
ter the last process, the number of 6.2 per cent sources in 
the sample-stars-reduced to 10492. The two-colour diagrams 
(m - g )o - {g - r )o and {g - r )o - (r - i )o for the 
final sample are given in Fig. 5. A few dozen stars with 
(w — 5 )o ~ —0.10 and {g — r )o ~ 0.20 are probably stars 
of spectral type A. 

4.3 Absolute magnitudes, distances, population 
types and density functions 

In the Sloan Digital Sky Survey (SDSS) photometry, the 
blue stars in the range 15 < g'* < 18 are dominated by thick- 
disc stars with a turn-off {g* — r*) ~ 0.33, and for g* > 18, 
the Galactic halo, which has a turn-off colour {g*—r*) ~ 0.2, 
becomes significant. Red stars ((?*— r*) > 1.3, arc dominated 
by thin-disc stars for all apparent magnitudes (Chen et al. 
2001). We used the same procedure to demonstrate the three 



populations (Fig. 3) and to determine the absolute magni- 
tudes for stars in each population by appropriate colour- 
magnitude diagrams. In our case, the apparent magnitude 
which separates the thick disc and halo stars seems to be a 
bit fainter relative to the SDSS photometry, i.e. g^ ~ 19, 
and the colour separating the red and bluer stars is slightly 
more blue, i.e. {g — r)o = l.l (Fig. 6). The absolute mag- 
nitudes of thick disc and halo stars are evaluated by means 
of the colour-magnitude diagrams of 47 Tuc ([Fe/H]=-0.65 
dex, Hesser et al. 1987) and M13 ([Fe/H]=-1.40 dex. Richer 
& Fahlman 1986) respectively, whereas for thin disc stars 
we used the colour-magnitude diagram of Lang (1992) for 
population I stars. The colours and absolute magnitudes in 
the UBV system were converted to ELAIS photometry as 
follows: We used two equations for colour transformations 
fi-om the WEB page of CASU^: 

g -r = 0.908(5 -V)+ 0.048, (12) 
t?' - B = -0.531(B - y) -h 0.053. (13) 

3 

http: // www. ast . cam. ac. uk / ~ wf csur /technical / photom / colours / 
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Figure 7. Spatial location for stars with absolute magnitude 6 < M{g ) < 7 as a function of apparent magnitude: (a)(17.0,17.5], (b) 
(17.5,18.0), (c)(18.0,18.5], (d)(18.5,19.0], (e)(19.0,19.5], (f)(19.5,20.0] and (g)(20.0,20.5]. The arrows correspond to the distances from the 
Galactic plane separating the populations (TN: thin disc, TK: thick disk and H: halo). 
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Figure 6. A [g — r')o colour histogram as a function of apparent 
magnitude, for the star sample: (a) for 17 < gg < 22, (b) for 
17 < < 19 and (c) for 19 < jiq < 22. The vertical downward 
arrow shows the limit value (g — r)o = l.l mag which separates 
the thin disc and the thick disc-halo couple. 



The first equation transforms B — V to g — r colour. Tlie 
second equation can be written in the following form which 
provides M{g ) absolute magnitudes: 



M{g ) = M{B) - 0.531(5 -V)+ 0.053. 



(14) 



The M{B) absolute magnitudes in eq. 14 were evaluated 
either by eq. 15 (for the data of Lang, 1992) or by eq. 16 
(for the data of clusters 47 Tuc and M13): 



M(B) = M{V) + (B-V), 
M{B) = B - {V - M{V))o, 



(15) 
(16) 



where {V — M{V))o is the distance modules of the cluster in 

question. The distance to a star relative to the Sun is carried 
out by the following formula: 



[g -M{g )]o = 51ogr-5. 



(17) 



The vertical distance to the Galactic plane (2) of a star 
could be evaluated by its distance r and its Galactic lati- 
tude (6) which could be provided by its right ascension and 
declination. 

The precise separation of stars into different popula- 
tions has been carried out by their spatial positions as a 
function of their absolute and apparent magnitudes. The 
procedure is based on the histograms for distance z from 
the Galactic plane. Fig. 7 gives the histograms for stars 
with & < M{g ) < 7 for the apparent g^ magnitude in- 
tervals 17 < g'o < 17.5, 17.5 < g'o < 18, 18 < po < 18.5, 
18.5 < 5o < 19, 19 < ffo < 19.5, 19.5 < go < 20 and 
20 < go ^ 20.5 as an example. The vertical arrows show the 
position that the number of stars decline. The distance z 
corresponding those positions are adopted as the borders of 
three populations of the Gala^xy, i.e. they limit the efficiency 
regions of the populations. 
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Figure 8. Absolute magnitude ranges dominated by different 
populations (panel a) and the break of these contributions down 
by distance bins for the absolute magnitudes 4 < M{g ) < 5, 5 < 
M(g') < 6, 6 < M{g') <7,7 < M{g') < 8 and 8 < M{g') < 13 
in panels (b), (c), (d), (e) and (f) respectively. Symbols: a plus 
denotes a halo, a square denotes a thick disc and a filled circle 
denotes a thin disc. 

This topic can be clarified in more detail as follows: 
The dominance regions of thin disc, thick disc and halo are 
at short, intermediate, and large z distances respectively. 
The number counts for thin disc decreases with increasing 
z distance, whereas for thick disc the number counts are 
very small at short z distances, then they increase at in- 
termediate z distances. The decreasing or increasing rate 
is not the same for the two populations. The same case is 
valid for thick disc and halo at large z distances. Hence, 
statistically, there should be drops in the number counts 
at population transitions. This technique first used in the 
paper of Karaali (1994). He showed that the ^-histograms 
for three populations show a multimodal distribution where 
the modes at short, intermediate and large z distances cor- 
respond to thin disc, thick disc and halo stars respectively. 



Table 2. Dominant regions for three populations: the thin disc, 

the thick disc and the halo, as a function of absolute and apparent 
magnitudes. The symbol z is the distance to the Galactic plane 
in kpc. 
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As shown by Karaali, the agreement of the kinematical dis- 
tribution of the sample stars with their spatial location is a 
strong confirmation of the technique in question. 

The technique improved in the recent years (cf. KBH) 
by introducing the apparent magnitude of stars used in the 
histograms. A population breaks at higher z-distances when 
one goes to faint apparent magnitudes, and there are his- 
tograms where the statistical fluctuations are rather small 
relative to the deeps at the population transitions. These 
two arguments are the clues in the separation of stars into 
different population types, i.e. thin and thick discs and halo. 
However, any wrong identification of a genuine drop reflects 
in the value of the parameter in question and its correspond- 
ing error. 

Table 2 gives the full set of absolute and apparent mag- 
nitude intervals and the efficiency regions of the populations. 
The distance over which a population, the thin disc for ex- 
ample, dominates increases with declining absolute magni- 
tude. That is, the three populations are not squeezed into 
small isolated volumes. The same holds also when one goes 
to apparently faint magnitudes in an absolute magnitude in- 
terval. These findings were cited also in our previous paper 
(KBH) and arc consistent with the results of Rcid & Ma- 
jewski (1993), who argued that the thick disc extends up to 
a ~ 4 kpc, a distance from the Galactic plane where halo 
stars cannot be omitted. Halo stars dominate the absolutely 
bright intervals, thick disc stars indicate the intermediate 
brightness intervals and thin disc stars indicate the faint in- 
tervals, as expected (Fig. 8a). If we break these contributions 
down by distance bins, we would reveal that the eflicient re- 
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Table 3. The logarithmic space density function D* = logD + 10, for different absolute magnitude intervals for the thin disc, r* = 
[(r'j + r2)/2]^/'^ is the cciitroid distance for the volume AVi_2. and z* = r* sin 6, (&) being the Galactic latitude of the field center. The 
other symbols are explained in the text (distances in kpc, volumes in pc"^). 
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Table 4. The logarithmic space density function, D* = log D + 10, for different absolute magnitude intervals for the thick disc. Symbols 
as in Table 3. 
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gion for each population shifts to shorter distances relative 
to the Sun, when one goes from absolutely bright to abso- 
lutely faint magnitudes (Fig. 8b-f). For example, thin disc 
is efficient at r ~ 1.5 kpc for the absolutely magnitudes 
6 < M{g') < 7 and 7 < M(g') < 8 whereas the efficiency 
shifts to r ~ 0.5 kpc for the interval 8 < M{g ) < 13. 

The logarithmic density functions, D* = log D + 10, are 
given in Tables .3-5 and Figures 9-11 for different absolute 
magnitudes for three populations, where: D = N/A.Vi,2; 
AVi,2 = (7r/180)^(n/3)(ri - r?); □ denotes the size of the 
field (6.571 dcg'^); ri and r2 denote the limiting distance of 
the volume AVi,2; A' denotes the number of stars (per unit 
absolute magnitude); r* — [{r\ + r|)/2]^''^ is the centroid 
distance for the volume AVi,2; and z* = r*sin6, (6) being 
the Galactic latitude of the field center. The horizontal thick 
lines, in Tables 3-5, corresponding to the limiting distance of 
completeness {zi) are evaluated by the following equations: 



Table 5. The logarithmic space density function, D* = lo 
10, for different absolute magnitude intervals for the halo, 
bols as in Table 3. 
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zi = n sin b, 



(19) 



[g -M{g )]o = 51ogn-5. 



(18) 



where (/() is the limiting apparent magnitude (17 and 20.5 
for the bright and faint stars respectively), n the limiting 
distance of completeness relative to the Sun and M{g ) is the 
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Figure 9. Comparison of the observed space density function with the density laws for different absolute magnitude intervals for the 
thin disc, (a) (5,6], (b) (6,7], (c) (7,8], (d) (8,9], (e) (9,10], (f) (10,11], (g) (11,12] and (h) (12,13]. The continuous curve represents the 
exponential law, the dashed curve, represents the sech law and dashed-dot curve represents the sech^ law. 




Figure 10. Comparison of the derived space density function 
with the exponential density law for different absolute magnitude 
intervals for the thick disc, (a) (4,5], (b) (5,6], (c) (6,7] and (d) 
(7,8). 



Figure 11. Comparison of the derived space density function 
with the de Vaucouleurs density law for different absolute mag- 
nitude intervals for the halo, (a) (4,5], (b) (5,6] and (c) (6,7]. 



appropriate absolute magnitude Mi or M2 for the absolute pic method. We have postulated mono-rnetallic stellar pop- 
magnitude interval Mi — M2 considered. ulations with no abundance gradient, and we have not ap- 
We acknowledge that in this work we have used a sim- plied any correction for binarism, contamination by compact 
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Table 6. Galactic model parameters for different absolute magnitude intervals for the thin disc resulting from the comparison of observed 
logarithiuic space densities with a (unique) density law (Fig. 9). The columns give: absolute magnitude interval M(g ), the density law, 
the logarithmic local space density n*, scaleheight for sech or sech^ density law H-^ (in the parenthesis), the scaleheight for exponential 
density law H, X^in' standard deviation s and the local space density for Hipparcos 0. 
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galaxies or giant /sub-giants neither. We are making a sim- 
plified evaluation. 



5 GALACTIC MODEL PARAMETERS 

Wc estimated Galactic model parameters by comparison of 
the derived space density functions, with the density laws 
both independently for each population as a function of ab- 
solute magnitude and simultaneously for all stellar popula- 
tions. 

5.1 Absolute magnitude dependent Galactic 
model parameters 

The thin disc density laws wore fitted with the additional 
constraint of producing local densities consistent with those 
derived from Hipparcos (Jahreiss k, Wielen 1997), a pro- 
cedure applied in our previous paper (KBH). It was dis- 
covered that the sech^ law fitted better for the intervals 
10 < M{g') < 11 and 11 < M{g') < 12 confirming our 
results in our paper mentioned above, however, contrary to 
our expectation, the exponential law fitted better for the 
interval 12 < M{g ) < 13 and the sech law fitted better 
for the interval 9 < M{g ) < 10, whereas the exponential 
law was favourite in our previous paper (KBH, Table 6 and 
Fig. 9). The comparison for absolutely bright intervals, i.e. 

8 < M(g') < 9, 7 < M{g') < 8, 6 < M{g') < 7 and 
5 < M{g )< 6 is carried out with the exponential law, as in 
our previous paper cited above. The scaleheight for thin disc 
increases monotonically from 163 to 363 pc when one goes 
from the absolute magnitude interval 12 < M{g ) < 13 to 
5 < M{g ) < 6, with the exception scaleheight H=211 pc for 
the interval 9 < M{g ) < 10 which is less than the one for 
the interval 10 < M(3 ) < 11, i.e. H=295 pc. As cited above, 

9 < M{g ) < 10 is the unique absolute magnitude interval 
where sech density law fitted better with the derived space 



densities in our work, and it is a transition interval between 
those for which either exponential law (for bright intervals) 
or sech^ law (for fainter intervals) fitted better. All scale- 
heights are equivalent to the exponential law scaleheights. 
The local space density for thin disc, for different absolute 
magnitude intervals, is consistent with the Hipparcos' one 
(Table 6). 

For the thick disc, the derived logarithmic space den- 
sity functions are compared with the exponential density 
law for the absolute magnitude intervals 7 < M{g ) < 8, 
6 < M{g') < 7, 5 < M{g') < 6 and 4 < M{g') < 5, (Ta- 
ble 7 and Fig. 10). The range for the scaleheight is rather 
small, 839-867 pc and the scaleheight itself is flat within the 
quoted uncertainties. The local space density relative to the 
local space density of thin disc (712 /ni) could not be given 
for the interval 4 < M{g ) < 5 due to lack of local space 
density for this interval for the thin disc. For the intervals 

5 < M{g') < 6 and 6 < M{g') < 7, n2/ni is 7.59 and 7.41 
per cent respectively, equivalent to the updated numerical 
values, whereas for the faintest interval, 7 < M{g ) < 8, 
712 /ni =3.31 per cent is close to the original value (Gilmore 

6 Wyse 1985). 

The derived logarithmic space density functions for 
the halo are compared with the de Vaucouleurs density 
law for the absolute magnitude intervals 6 < M(g ) < 7, 
5 < M{g') < 6 and 4 < M{g') < 5 (Table 8 and Fig. 11). 
The local space density relative to the thin disc (ns/ni) 
could not be given for the interval 4 < M{g ) < 5 due to 
the reason cited above. For the intervals 6 < M(g ) < 7 
and 5 < M{g ) < 6, n3/ni=0.04 and 0.06 per cent respec- 
tively. The numerical values for the axial ratio k for two 
intervals are close to each other, i.e. k = 0.78 and k — 0.73 
for 6 < M{g') < 7 and 5 < M{g') < 6 respectively, but a 
bit less for the interval 4 < M{g ) < 5, k = 0.60, consistent 
with the previous ones within the uncertainties however. 

The parameters derived for three populations have been 
tested by the luminosity function (Fig. 12), where ip*{M) 
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Table 7. Galactic model parameters for the thick disc. n2/ni 
indicates the local space density for the thick disc relative to the 
thin disc. Other symbols are same as in Table 6. 
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Table 8. Galactic model parameters for the halo, k and ns/ni 
give the axial ratio and the local space density for the halo relative 
to the thin disc, respectively. Other symbols axe as in Table 6. 
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10.21 
0.14 
21.08 


± 0.26 
0.01 
0.24 


0.04 
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is the total of the local space densities for three popula- 
tions. The local space densities for the thick disc and the 
halo arc presented in Table 9. The local space densities of 
Hipparcos were converted to ELAIS colours by the combi- 
nation of eqs (13) and (15) which give the following relation 
between M{g ) and M{V) absolute magnitudes: 

M{g') = M{V) + 0.469(5 -V)+ 0.053 (20) 

There is a good agreement between our luminosity func- 
tion and that of Hipparcos ( Jahreiss & Wielen 1997) . Also 
the error bars are rather short. We used the procedure of 
Phleps et al. (2000) for the error estimation in Tables 6-8 
(above) and Tables 13 and 14 (in the following sections), i.e. 
changing the values of the parameters until Xmin increases 
or decreases by 1. 

5.2 Model parameter estimation by simultaneous 
compcirison to the Galactic stellar populations 

We estimated the model parameters for three populations si- 
multaneously by comparison of the combined derived space 
density functions with the combined density laws. We car- 
ried out this work for two sets of absolute magnitude inter- 
vals, 4 < M{g') < 10 and 4 < M{g') < 13. The fit for the 
second interval was done due to our experience that model 



Table 9. Local luminosity functions for thick disc {ip''(M)xk) 
and halo {(p*{M)[j). The luminosity function oi Hipparcos is also 
given in the last column. 
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4.19 


7.47 
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5.66 


4.43 
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Figure 12. The local luminosity function obtained from com- 
bining the local space densities for the thin and thick discs and 
the halo, resulting from comparison of the derived space density 
function with the density laws, for different absolute magnitude 
intervals. The "0" symbols show the Hipparcos values. TK and 
H corresponds to only thick disc and halo local luminosity, re- 
spectively. 

parameters are absolute magnitude dependent and that it 
covers the thin-disc stars with 10 < M{g ) < 13, the density 
functions of which behave differently from the density func- 
tions for stars with other absolute magnitudes. Actually, we 
will see in the following that the luminosity function for the 
absolute magnitude interval 4 < M{g ) < 13 differs from 
the one for stars with 4 < M{g ) < 10 considerably (Fig. 14 
and Fig. 16) . The number of stars as a function of distance r 
relative to the Sun for nine absolute magnitude intervals are 
given in Table 10 and the density functions per unit abso- 
lute magnitude interval evaluated by these data are shown 
in Tables 11 and 12 for the intervals 4 < M{g') < 10 and 
4 < M{g ) < 13 respectively. 

5.2.1 Model parameters by means of absolute magnitudes 
4 < M{g') < 10 

The combined derived densities per absolute magnitude in- 
terval for three populations, the thin and thick discs and 
the halo for stars with 4 < M{g ) < 10 (Table 11), are com- 
pared with the combined density laws (Fig. 13). The derived 
parameters are given in Table 13. All these parameters are 
in agreement with the ones given in Table 1, and they lie 
between two corresponding parameters cited in Section 5.1, 
except the scaleheight of thick disc, 760 pc, which is rather 
smaller than the scaleheight 839 pc, the smallest one in Table 
7. Thus, the scaleheight 269 pc for the thin disc is between 
the ones for absolute magnitude intervals 10 < M{g ) < 11 
and 11 < M{g ) < 12, and the logarithmic local space den- 
sity 7.51 is almost equal to the corresponding one for the 
absolute magnitude interval 9 < M{g ) < 10. For the thick 
disc, the local space density relative to thin disc, 6.46 per 
cent, is between the local space densities for the absolute 
magnitude intervals 6 < M{g') < 7 and 7 < M{g') < 8. It 
is interesting that the local space density for the halo relative 
to thin disc, 0.08 per cent, is almost equal to the one for the 
absolute magnitude interval 5 < M{g ) < 6 in Table 8, how- 
ever, the axial ratio k = 0.55 is considerably smaller than 
the ones in the same table. The resulting luminosity func- 
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Table 10. Number of stars as a function of distance r relative to the Sun for nine absolute magnitude intervals (distances in kpc). 
Horizontal thick lines correspond the limiting distance of completeness. 
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Figure 13. Comparison of the derived and combined space den- 
sity function for the thin and thick discs and the halo with the 
combined density law, for stars with 4 < M(g ) < 10. 



tion (Fig. 14) from the comparison of the model with these 
parameters and the combined derived density functions per 
absolute magnitude interval is in agreement with the one 
of Hipparcos (Jahreiss & Wielen 1997). However, the error 
bars are longer than the ones in Fig. 12, particularly for the 
faint magnitudes. 




mi 

Figure 14. The local luminosity function resulting from the com- 
parison of the combined derived space density function with the 
combined density law, for stars with 4 < M{g ) < 10. 



Sun relative to stars brighter than M{g ) = 10, and they 
affect the combined density function considerably. Also the 
corresponding luminosity function is not in agreement with 
the one of Hipparcos (Fig. 16). 



5.2.2 Model parameters by means of absolute magnitudes 
4 < M{g') < 13 

We carried out the work cited in previous paragraph for 
stars with a larger range of absolute magnitude, i.e. 4 < 
M{g ) < 13. The derived density function is given in Ta- 
ble 12 and its comparison with the combined density law 
is shown in Fig. 15. Most of the derived parameters (Ta- 
ble 14), especially the local densities, are rather different 
than the ones cited in Section 5.1 and 5.2.1. The reason 
for this discrepancy is that stars with absolute magnitudes 
10 < M{g ) < 13 have relatively larger local space densities 
{Hipparcos; Jahreiss & Wielen 1997) and are closer to the 



6 DISCUSSION 

We estimated the Galactic model parameters by comparison 
of the derived space density functions per absolute magni- 
tude interval, in the perpendicular direction to the Galactic 
plane, with a unique density law for each population individ- 
ually for the ELAIS field (a = 16''10'"00'', d = +54°30'00"; 
/ = 84°. 27, b = +U°.90; 6.571 deg^ epoch 2000), by Vega 
photometry. The separation of stars into different popular- 
tions has been carried out by their spatial position as a 
function of both absolute and apparent magnitude (KBH, 
see also Karaali 1994). This work covers nine absolute mag- 
nitude intervals, i.e. 4 < M{g') < 5, 5 < M{g') < 6, 
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Table 11. The logarithmic space density function, D* = logD + 
10, per unit absolute magnitude interval for stars with 4 < 
M{g ) < 10. < W > is the mean of number of stars weighted 
by their numbers. The data are taken from Table 10. The other 
symbols are explained in the text (distances in kpc, volumes in 

pc3). 
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Table 12. The logarithmic space density function, D* = logD-|- 
10, per unit absolute magnitude interval for stars with 4 < 
^{g') < 13 (symbols as in Table 11). 



n — r2 


AVi,2 


r* 


2* 


<N> 


D* 


0.2-0.4 


3.74 


(4) 


0.33 


0.23 


152 


7.61 


0.4-0.7 


1.86 


(5) 


0.59 


0.42 


293 


7.20 


0.7-1.0 


4.38 


(5) 


0.88 


0.62 


221 


6.70 


1.0-1.5 


1.58 


(6) 


1.30 


0.92 


230 


6.16 


1.5-2.0 


3.09 


(6) 


1.78 


1.26 


223 


5.86 


2.0-2.5 


5.09 


(6) 


2.28 


1.61 


180 


5.55 


2.5-3.0 


7.59 


(6) 


2.77 


1.96 


170 


5.35 


3.0-4.0 


2.47 


(7) 


3.57 


2.52 


277 


5.05 


4.0-5.0 


4.07 


(7) 


4.56 


3.22 


213 


4.72 


5.0-7.5 


1.98 


(8) 


6.49 


4.58 


212 


4.03 


7.5-10.0 


3.86 


(8) 


8.92 


6.30 


140 


3.56 


10.0-12.5 


6.36 


(8) 


11.39 


8.04 


90 


3.15 


12.5-15.0 


9.49 


(8) 


13.86 


9.78 


85 


2.95 


15.0-17.5 


1.32 


(9) 


16.35 


11.54 


41 


2.49 



6 < M{g') < 7, 7 < M{g') < 8, 8 < M{g) < 9, 
9 < M(g ) < 10, 10 < M(g') < 11, 11 < M{g') < 12 
and 12 < M{g ) < 13. However, the populations are not 
dominant in all absolute magnitude intervals. We consider 
two density laws for the thin-disc stars: the density law 
sech^ fits better with the derived space densitp^ functions 
for the absolute magnitude intervals 10 < M{g ) < 11 and 
11 < M{g ) < 12, whereas for the absolute magnitude in- 
tervals 5 < M{g') < 6, 6 < M{g') < 7, 7 < M{g) < 8 and 
8 < M{g ) < 9 the exponential density law is favorable as in 
our previous paper (KBH). Contrary to our expectation, for 
the absolutely faintest interval 12 < M(p ) < 13 the expo- 
nential law (not the sech^ one) fits better with the derived 
space density functions. 9 < M{g ) < 10 is a transition ab- 
solute magnitude interval where sech fits better, a density 
law which does not correspond to anything physical. The 
scaleheight for the thin disc increases monotonically from 
163 to 363 pc, when one goes from the absolute magnitude 



Table 13. Galactic model parameters estimated by comparison of 

the logarithmic space density function for stars with 4 < M(g ) < 
10 (given in Table 11) and the combined density laws (Fig. 13). 
Symbols; n* is the logarithmic space density, H is the scaleheight, 
n/ni is the local space density relative to the thin disc and k is 
the axial ratio for the halo. 
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Table 14. Galactic model parameters estimated by comparison of 
the logarithmic space density function for stars with 4 < Af (g ) < 
13 (given in Table 12) and the combined density laws (Fig. 14). 
Symbols are as in Table 13. 
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interval 12 < M(g') < 13 to 5 < M{g') < 6 with the 
exception scaleheight H=211 pc for the transition interval 

9 < M{g ) < 10 which is less than the one for the interval 

10 < M{g ) < 11, i.e. H=295 pc. Some researchers restrict 
their works related with the Galactic model estimation to 
absolute magnitude. The recent work of Robin ot al. (2003) 
who treated stars with M{V) < 8 can be given as an exam- 
ple. Thus, if we consider the range of the scaleheight only 
for stars with 4 < M(g' ) < 10 in our work, we notice that 
it almost overlaps with the range of the scaleheight defined 
by the minimum and maximum scalclioiglits given in Table 
1, i.e. 200-350 pc (we did not take into account the upper 
limit for the thin disc of Robin & Creze (1986), 475 pc). The 
local space density for the thin disc decreases monotonically 




5 10 15 

z(kpc) 



Figure 15. Comparison of the derived and combined space den- 
sity function for the thin and thick discs and the halo with the 
combined density laws, for stars with 4 < M(p ) < 13. 



© 0000 RAS, MNRAS 000, 1-?? 



14 Bilir et al. 




cited in the literature (see Table 1) and the larger number 
of stars results the larger local space density. Hence, we can 
add another point; and 

(iii) The largeness of the local space density for an abso- 
lute magnitude interval for a specific population is propor- 
tional to the dominance of that population in the interval 
considered. Therefore, if a population is dominant in an ab- 
solute magnitude interval, then the local space density for 
that population is large. 
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Figure 16. The local luminosity function resulting from the com- 
parison of tiic combined derived space density function with the 
combined density law, for stars with 4 < M{g ) < 13. 



within the uncertainties, from absolutely faint magnitude 
intervals to the bright ones, however, the gradient converges 

to zero at the bright intervals in agreement with the cor- 
responding local densities of Hipparcos (Jahreiss & Wielen 
1997). 

The logarithmic space densities for the thick disc could 
be derived for the absolute magnitude intervals 4 < M(p' ) < 
5, 5 < M(g') < 6, 6 < M(g') < 7 and 7 < M{g') < 8. The 
range for the scaleheight is rather small, 839-867 pc, and the 
scaleheight itself is flat within the quoted uncertainties. The 
scaleheight for the thick disc in our work is in good agree- 
ment with the scaleheight claimed by many authors (cf. Ya- 
magata & Yoshii 1992, von Hippcl & Botlmn 1993, Buscr ct 
al. 1998; 1999). However, contrary to our expectation, the 
range of the scaleheight for the thick disc cited very recently 
is large: actually Chen ct al. (2001) and Siegel et al. (2002) 
give 0.58-0.75 and 0.7-1.0 kpc respectively. Now, let us com- 
pare these ranges with the ones in Tables 1.3 and 14, where 
the Galactic model parameters estimated by simultaneous 
comparison to the Galactic stellar populations correspond 
to stars with 4 < M{g') < 10 and 4 < M{g') < 13 respec- 
tively. The range for the scaleheight of thick disc is 705-822 
for Table 13 and 737-921 pc for Table 14. We can easily see 
that they are in good agreement, especially if we round our 
values, i.e. 0.7-0.9 kpc, they overlap with the ones of Siegel 
ct al. (2002). We should remind that the absolute magni- 
tudes of stars treated by Siegel et al. (2002) extend down 
to M{R) = 10.2 which corresponds to M{g') > 10.2 mag 
in the Vega photometry. This comparison encourage us to 
argue two points: 

(i) The absolutely faint stars cause large ranges in the 
estimation of Galactic model parameters; 

(ii) Galactic model parameters are mass (and hence abso- 
lute magnitude) dependent as claimed in our previous paper 
(KBH). The local space density for the thick disc relative 
to thin disc increases from 3.31 to 7.59 per cent, from abso- 
lutely faint magnitude intervals to the bright ones. The large 
values correspond to absolute magnitude intervals where 
thick disc is dominant (Table 7, Fig. 8). Here we reveal an- 
other property: although the range of absolute magnitude 
for the thick disc stars, 4 < M{g ) < 8, is not large, the lo- 
cal space density for the thick disc cover almost all the range 



The logaritmic space density functions for the halo 
could be derived only for throe absolute magnitude inter- 
vals, 4 < M{g') < 5, 5 < M{g') < 6 and 6 < M{g') < 7, 
and they are compared with the de Vauculeurs density law. 
The local space density relative to the thin disc (na/ni) 
could not be given for the interval 4 < M{g ) < 5, due to 
the lack of local space density for this interval for the thin 
disc. For the intervals 5 < M{g') < 6 and 6 < M{g') < 7, 
ns/ni = 0.06 and ns/ni = 0.04 respectively, consistent with 
the results of Buser ct al. (1998, 1999) and KBH. The rm- 
merical values for the ax;ial ratio k for two intervals are close 
to each other, i.e. k = 0.73 and k = 0.78 for 5 < M{g') < 6 
and 6 < M{g ) < 7 respectively, but a bit less for the inter- 
val 4 < M{g ) < 5, K = 0.60, consistent with the previous 
ones within the uncertainties however. It is interesting that 
the axial ratio estimated by simultaneous comparison to the 
Galactic stellar populations for stars within 4 < M{g ) < 13 
is lesser, k = 0.55 (Tables 13 and 14) than the ones cited 
above. It seems that absolutely faint magnitudes where halo 
stars are very rare causes bad k estimation. 

Finally, we compared the Galactic model parameters 
estimated in our previous work (Z = 68°. 15, b = —48°. 38) 
and in this work {I = 84°. 27, b = -|-44°.90) in Table 15. The 
results confirm the idea that the Galactic model parameters 
are mass (and hence absolute magnitude) dependent. There 
is a good agreement between two sets of data, however, we 
should note some points and keep in mind in the comparison 
with the results that would appear in future: 

(a) Although the scaleheights for a specific abso- 
lute magnitude interval of thin disc for two fields are rather 
close to each other, the local space density for the inter- 
val 11 < M{g') < 12 is a bit larger for the field SA 114 
{nl = 8.6) than the one for the ELAIS field {nl = 8.0). 
This slight discrepancy probably originates from the previ- 
ous work because the corresponding total local space density 
does not agree with the Hipparcos (Jahreiss & Wielen 1997) 
one either; 

(b) Although the scaleheights estimated for the thick 
disc for two fields are rather close to each other, the local 
space density relative to the local space density of thin disc 
for the field SA 114 is larger than the one for ELAIS field; 
and 

(c) The axial ratio, k, for the halo for two fields are 
almost the same, whereas the local space density relative the 
local space density of the thin disc for the absolute magni- 
tude interval 5 < M{g ) < 6 for the field SA 114 is 2.5 times 
that of the corresponding one for ELAIS field. 
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Table 15. Comparison of the most appropriate Galactic model parameters for two works: (SA 114) for our previous work (KBH) and 
(ELAIS) for the present work. 
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